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ABSTRACT 

We report here X-ray imaging spectroscopy observations of the northeastern shell of the supernova remnant 
RCW 86 with Chandra and XMM-Newton. Along this part of the shell the dominant X-ray radiation mechanism 
changes from thermal to synchrotron emission. We argue that both the presence of X-ray synchrotron radiation 
and the width of the synchrotron emitting region suggest a locally higher shock velocity of V s ~ 2700 km s _1 
and a magnetic field of B rs 24 ± 5 fiG. Moreover, we also show that a simple power law cosmic ray electron 
spectrum with an exponential cut-off cannot explain the broad band synchrotron emission. Instead a concave 
electron spectrum is needed, as predicted by non-linear shock acceleration models. Finally, we show that the 
derived shock velocity strengthens the case that RCW 86 is the remnant of SN 185. 
Subject headings: shock waves - X-rays: observations individual (RCW 86) - supernova remnants 
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1. INTRODUCTION 

Since t he discovery of X-ra y synchrotron radiation from 
SN1006 i lKovamaetalJll995l) it has been found that most 
young, shell-type sup ernova remn ants (SNRs) emit X-ray 
synchrotron radiation i Vinkl 1200a. for a review). For the 
youngest SNRs Cas A, Kepler (SN1604) and Tycho (SN1572) 
this radiation is confined to a narrow region close to the shock 
front. This does not seem to be the case for somewhat older, 
but physically m uch larger objects lik e RCW 86, G 266.2- 
1.2 Jsiane et alJl2001l) and G347. 3-0.5 (ICassam-Chenai et alJ 
120041) 7 The X-ray emission from the latter two SNRs does 
in fact only reveal synchrotron rad iation, and both have be en 
detected in TeV gamma-rays ( Aharonian et al l27Ml200l . 

It has been shown that the thickness of the X-ray syn- 
chrotron emitting region is directly related to the post-shock 
magnetic field strength (e.g lVink & Lam ing 2003). The rel- 
atively strong fields found are probably a result of m agnetic 
field amplification by co smic ray streaming (e.g. lBellll2004t 
iBvkov & Toptvgirl20 05). The presence of X-ray synchrotron 
radiation is in itself an indication of efficient acceleration, 
and if the cut-off energy of the photon spectrum is deter- 
mined by synchrotron losses, it is independent of the mag- 
netic field strength, but scal es with shock velocity, as oc 
(Aharonian & Atov anl 19991) . 

Here we report on the analysis of Chandra and XMM- 
Newton data of the northeastern (NE) pail of RCW 86 
(G315.4-2.1). RCW 86 is an interesting SNR. The non- 
thermal X-ray emitting regions are broader than those of the 
historical SNRs, and not confined to the for ward shock region , 
possibly as a result of projection effects SV'mk et alJll997l) . 
In that respect RCW 86 resembles G266.2-1.2 and G347.3- 
0.5. However, unlike those SNRs RCW 86 also emits notice- 
able thermal X-ray emission. This allows the determination 
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of plasma properties, which can help to determine the condi- 
tions that may potentially lead to X-ray synchrotron emission 
in older SNRs. In particular the NE part seems best suited 
for that purpose, since the X-ray synchrotron radiation is con- 
fined to the region directly behind the shock front, and, as a 
result, the geometry of the emitting region is easier to asses. 

2. OBSERVATIONS AND ANALYSIS 

The NE of RCW 86 was observed by Chandra and XMM- 
Newton for their joint program. The Chandra observation 
(ID 500436) has a net exposure time of 72.6 ks and was 
made on June 6, 2004. The XMM-Newton observation (ID 
0208000101) was made on January 26, 2004 with an expo- 
sure time of 60 ks. Both Chandra and XMM- Newton ob - 
served other parts of RCW 86 before fe.g. iRho et afll2002ft . 
Here we present Chandra imaging of the NE part of RCW 86 
(Fig.|4]Plate 1), but for spectrosco py we focus on the spectra 
obtained with the EPIC-MOS (iTurner et alJl2001l) instrument 
of XMM-Newton, since it offers the best spectral resolution. 
For comparison we also present an analysis of EPIC-MOS 
data of an observation of the northern part of RCW 86 (Obs 
ID 0110011401, 18 ks). Moreover, we used all six observa- 
tions of RCW 86 by XMM-Newton to produce the X-ray map 
shown in Fig.^ It also shows the spectral extraction regions. 
Background spectra were taken from regions outside the SNR. 
The data reduction was done with the standard software pack- 
ages ciao 3.0 for Chandra and SAS 6.5.0 for XMM-Newton. 

The spectrum from the region labeled "NE" (Fig. 
shows an absence of bright line feature s. This is u sual! 
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taken as a sign for synchrotron emission (Bamba et al. 200( 
iBorkowski etaDl2001l IRho et alJ l2002). but for RCW 86 it 
has been argued by some of us in the past that non-thermal 
brem sstrahlung may be viable mechanism as well dVink et all 
fl997h . One of the reasons was the detection of Fe-K emis- 
sion at 6.4 keV in the southwest coinciding with t he hard X- 
ray co ntinuum. This interpretation was disputed bv lRho et alJ 
(2002), based on energetic grounds. Here we concur with that 
view. In fact the XMM-Newton and Chandra spectra of the NE 
region do not show any evidence for Fe-K emission, providing 
additional evidence for an X-ray synchrotron interpretation. 

Figure |4] shows that along the shell the X-ray emission 
changes from thermal radiation to predominantly synchrotron 
radiation. Unlike in other parts of RCW 86, it is unambiguous 
that in the NE the non-thermal emission starts at the shock 
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FIG. 1. — Left: XMM-Newton-EPIC (MOS/PN) mosaic of RCW 86, with a color coding similar to Fig.|4](Plate 1), using ajsquare root brightness scaling 
Spectral extraction regions are overlayed. Right: archival MOST 0.84 GHz radio map I Whiteoak & Green 1996; Dickel et al. 2001). 



TABLE 1 

Best fit models FOR THE XMM-Newton MOS 1 &2 SPECTRA. 





NE 


E bright E faint 


N bright 


N faint 


EMi (10 53 cm" 3 kpc- 2 ) 


12.6 ± 1.9 


311 ±45 70.6 ±8.0 


219 ± 48 


69.5 


kT cl (keV) 


6.7±2.6 


0.57 ±0.05 0.93 ±0.07 


0.96 ±0.13 


6.3 ± 1.6 


n c ti (10 9 cm" 3 s) 


2.25 ±0.15 


6.7 ±0.6 5.67 ±0.03 


4.0 ±0.3 


2.27 ±0.12 


EM 2 (10 53 cm" 3 kpc[- 2 ) 




43.8 ±5.3 


40.9 ± 5.6 




fcT e2 (keV) 




3.0 ±0.3 


3.2 ±0.6 




n a t 2 (10 9 cm- 3 s) 




17.0 ±0.5 


19.7 ±0.9 




O 




0.56 ± 0.06 


0.68 ±0.09 




Ne 




0.59 ± 0.04 


0.72 ±0.07 




Mg 




0.44 ± 0.04 


0.39 ±0.05 




Si 




0.27 ±0.04 


0.38 ±0.07 




Fe 




0.64 ±0.07 


0.79 ±0.14 




PLNorm(10- 3 s- 1 keV- 1 cm" 2 ) 


1.22 ±0.03 






1.8 ± 0.1 


r 


2.82 ±0.04 






3.5 ± 0.2 


Nu (10 21 cm" 2 ) 


4.1 ±0.1 


3.6 ±0.2 


3.6 ±0.2 


3.9 ± 0.2 


C-stat/d.o.f. 


190.6/169 


(302+260)/230 


182/85 


242/103 



NOTE. — The models consist of eit her one or two NEI compo nents, or of an NEI component with an additional power law component. Abudances are given 
with respect to the solar abunances of Anders & Grevesse 1 1989). An entry "1" means the abundance as fixed to the solar value. Entries covering two columns 
were obtained by fitting two spectra simultaneously, and forcing the parameters to be identical. The emission measure, EM, is defined as J n c nndV/ d 2 . The 
power law normalization is given at 1 keV. All errors are la errors (Ax 2 = 1). 



front. Surprisingly, the radio map shows that in the NE the 
X-ray synchrotron radiation is not associated with bright ra- 
dio emission, but instead relatively weak radio emission. In 
order to better determine the plasma conditions of the X-ray 
synchrotron emitting regions, we do not present here the spec- 
trum extracted from the whole X-ray synchrotron emitting re- 
gion, but from a smaller region, labeled "NE" in Fig.^ which 
shows more traces of thermal emission. Apart from the NE 
shell, we also analyzed a faint region outside the main shell 
("N-faint"), which shows up better in radio than in X-ray. 

We fitted the EPIC -MOS spectra using the spectral anal- 
ysis software SPEX ( iKaastra & Mewell2006T) . The spectra 



are reasonably well described by a combination of two non- 
equilibrium ionization (NEI) models or, for the X-ray syn- 
chrotron dominated regions, by one NEI model and a power 
law spectrum (Table The spectral fits indicate extreme 
departures from collisional equilibratio n, even compar ed to 
other young SNRs, except for SN1006 ( Vink et alJ2003l) . Al- 
though the fit quality is not always satisfactory, both the fit- 
ted models and the raw spectra itself (Fig. |2ji show that the 
fainter the thermal emission, the lower the ionization degree. 
This can be judged from the increase of O VIII with respect 
to the O VII emission and presence of Fe L emission from the 
brighter parts. There is also a hint that the X-ray synchrotron 
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FIG. 2. — XMM-Newton EPIC-MOS spectra from the regions labeled in Fig.Q Left: logarithmic plots of a thermal and a non-thermal spectrum. Right: a 
comparison of the line emission from various regions. From the northeastern spectrum (in red, left panel) the best fit power law model has been subtracted in 
order to emphasize the thermal emission. Dashed lines indicate (from left to right) the energies of OVII Hea, OVIIILya, and Fe XVIII line emission. 
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FIG. 3. — Broad band (uF v ) spectrum of the NE region with the 847 MHz 
flux density based on archival MOST data. The models consist of synchrotron 
radiation from power law, and concavely shaped electron energy distributions 
with exponential cut-offs at high energies, c onvolved with the synchrotron 
emission function i Ginzbure & Svrovatski 1969). The radio spectral index is 
fixed to a = —0.6 ICaswell et al!1975l) 



emitting regions contain hotter plasma, and are characterized 
by lower ionization parameters (n e t). However, we caution 
that, since most of the continuum is non-thermal, the tem- 
perature is solely determined by line ratios, making it more 
sensitive to potential errors in the atomic data. This may also 
explain the peculiarity that whenever a second thermal com- 
ponent was fitted, its best n c t value was much higher than 
the principal component. The second component improved, 
in particular, the fit between 1 and 1.2 keV, dominated by 
n = 3 — > 1 Ne IX line emission. 

3. DISCUSSION 

A recent development in the interpretation of X-ray syn- 
chrotron radiation is that the width of the region can be used 
to estimate the downstream magnetic field. However, dif- 
ferent groups have used, at fac e value, different meth ods 
for estimating the magnetic field: IVink & Lamind J2003I) as- 



sumed that the width, I, of the X-ray synchrotron emitting 
region is determined by a combination of plasma velocity rel- 
atively to the shock front, u, and the synchrotron loss time 
floss — 637/ B 2 E , i.e. l a d v — uti oss , with E the particle 
energy in erg and B the magnetic field strength in Gauss. The 
other method ass umes that the width corresponds to th e diffu- 
sion length scale (Bamb alTaTl27)0llVolk et all200fl . given 
by Idiff — D/u, with D = cE/3eB the diffusion coeffi- 
cient. In essence Idiff is the length scale at which advection 
starts to dominate over diffusion as means of transporting par- 
ticles. Both methods give a combination of E and B, which 
can be solved by using the fact that the observed photon ener- 
gies peak around e = 7 ABE 2 kcV.Both methods rely on dif- 
ferent assumptions. First of all, for standard shocks u ~V S /A, 
but for very efficient shock acceleration the compression fac- 
tor may be stronger than a factor 4. Secondly, the diffusion 
length method assumes D = cE/3eB, which is the smallest 
diffusion coefficient possible, as the particle mean free path is 
then equal to the gyro-radius (the "Bohm limit"). 

It turns out that with both metho ds very similar mag netic 
field estim ates are ob tained jBalletl 120051 I Vinkl 120061) . As 
shown bv lVinkl i2006l) . this is to be expected if one observes 
the X-ray synchrotron spectrum near the spectral cut-off ener- 
gies: The acceleration time for particles according to the first 
order Fermi accele ration theory is, within a factor of order 
one, r acc « D/u 2 (Malkov & Drurv 20011). For electrons the 
acceleration is limited by synchrotron losses. So there is only 
a net acceleration if r QCC < ti„ rs , an d the maximum energy is 
reached when r QCC « ti oss (Reynolds 1998). So: 

Tacc ~ Tioss D / 'u S3 UTi oss l dl ff « l a dv (1) 

Therefore, Idiff ~ ladv is the geometrical equivalent of 
Tioss ~ T acc . Note that the observational fact that the dif- 
fusion length scale and advection length scale methods give 
similar results is a justification for the assumption that the 
diffusion coefficient is close to the Bohm li mit, and a co m- 
pression ratio close to the standard value of 4 ( Vink 2006). 

The Chandra image reveals a width of the X-ray syn- 
chrotron shell of ~ 100", c orresponding to 3 . 7 x 10 18 cm 
for a distance of d = 2.5 kpc ( Wes terlundll969tlRosado et alJ 
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H99l . Fitting a projected shell model we estimate a phys- 
ical width of 1.7 x 10 18 cm. In addition we apply a fac- 
tor 0.6, because the actual width is a convolution of advec- 
tion and diffusion processes (i.e. we set Idiff = Udv in Eq. 
(1) of iBerezhko & Volkl J2004l) l. so l adv = 1.0 x 10 18 cm. 
Assuming that the sho ck velocity is V s = 600 kms" 1 
JGhavamian et all2001l) gives inconsistent results for the two 
methods: B ~ 90 /iG employing the diffusion length method, 
and 6 /iG assuming the advection length method. 

However, if we no longer assume a shock velocity of 
600 kms -1 , which after all was based on optical observa- 
tions of a different part of RCW 86 from which no X-ray syn- 
chrotron radiation is emitted, we can use the condition that 

ladv — Idi 

//in order to estimate V s and B, using an observed 
photon energy of e ks 1 keV. We found that B can be directly 
estimated from the diffusion/advection length alone, whereas 
the plasma velocity only depends on the photon energy: 



B 



, c , 
-Ye' 



1/3,-2/3 



adv 



= 24 



I, 



dv 



-2/3 



1.0 • 10 18 cm y 



MG,(2) 



V. 



— = 2650-£ 

7.4 3e • 637 4 V 1 keV 



km/s,(3) 



where \ the shock compression factor. E q. Q was previously 
reported bv lAharonian & Atovanl (11999ft . To get more accu- 
rate estimations one has to solve the, model dependent, ki- 
netic equations for the particle distribution. For the moment, 
we estimate the error in the physical width due to distance un- 
certainties and projection effects to be ~ 30%, resulting in an 
error in B of ~ 5 /iG. The downstream magnetic field there- 
fore appears to be lower than for other SNRs (IVinkl 2006). 
There is some uncertainty in the actual cut-off energy of the 
synchrotron spectrum (see below), but uncertainties about the 
assumptions - larger compression ratio, less efficient diffusion 
- make Eq. Q effectively a lower limit. 

One may wonder why the shock velocity in some regions 
may be so much higher, and why it has not resulted in a more 
distorted shell. Here the thermal spectra helps to answer this 
question. For all fitted regions the n c t value of the primary 
component is low, being n c t~ 6.7 x 10 9 cm~ 3 s even for the 
brightest region. From the size of our extraction box, we es- 
timate that the emitting volume is V = 10 56 — 10 57 cm 3 , 
together with the emission measure this implies an electron 
density of n c = 0.5 — 1.6 cm~ 3 . If we use this to estimate 
how long ago the plasma was shock heated we find t < 425 yr. 
This is surprisingly short for a large SNR as RCW 86. In 
such a time the difference in radius between regions with low 
and high shock velocity would be ~ 0.9 pc on an average 
shock radius of 16 pc. The short interaction time supports 
the idea that RCW 86 is a SNR expanding in a wind blown 
bubble jVinketalJll997h . Such SNRs expand rapidly for a 



long time, but the shock velocity drops rather suddenly as 
soon as the shock starts interacting with the surr ounding shell 
swept up by the stellar wind (e.g. Dwarkadas 2005). This 
suggests that in RCW 86 the shock has reached in some re- 
gions the dense shell around the bubble some 400 yr ago, af- 
ter which it rapidly decelerated. In other regions the shock 
velocity is still high, but, due to the low density, the ther- 
mal X-ray emission is weak. This explains the coexistence 
of relatively weak radio synchrotron emission with conspic- 
uous X-ray synchrotron emission: due to the lower density 
fewer electrons are accelerated, but because of the high shock 
velocity they can be accelerated to higher energies. 

Interestingly, the X-ray synchrotron radiation is relatively 
bright, with respect to the radio emission, because we find 
that the simplest broad band synchrotron model, i.e. syn- 
chrotron radiation from a power law electron spectrum with 
an exponential cut-off does not fit the data (Fig. [3}. It can 
explain the X-ray flux, but not the spectral slope. Instead the 
electron spectrum needs to be concave, as predicted by non- 
linear shock acceleration models. However, we can not deter- 
mine whether the electron power law index bends toward -2, 
predicted for an overall shock compression ratio of 4, or to- 
ward -1.5, predicted for str ongly cosmic ray modified shocks 
jBerezhko & Ellison! 19991) . 

The NE region of RCW 86 has properties resembling those 
of the TeV emitting SNRs G347. 3-0.5, and G266.2-1.2: weak 
radio emission, and X-ray emission (almost) entirely consist- 
ing of synchrotron radiation. For RCW 86 the (weak) ther- 
mal X-ray emission indicates that these properties are due 
to a low density combined with a, relatively, high shock ve- 
locity. We speculate that for G347.3-0.5 and G266.2-1.2 the 
shock also moves through a low density region, e.g. a stel- 
lar wind bubble, and the shock velocity is similarly high. 
RCW 86 may be different in that some parts of the shock 
have reached the shell. Inside the bubble the shock evolu- 
tion can be approximated by the Sedov self-similar model, 
but this breaks down as soon as the shell is reached. Finally, 
assuming a Sedov evolution and using the apparent radius of 
22', the age of the remnant is estimated to be t = 7^7- » 

2250(14/2700 km s^ 1 ) yr. This would put the explosion date 
of RCW 86 closer to A D 185, the year a put ative supernova 
was observed in China (Stephenson & Green 2002). A shock 
velocity of ~ 600 kms -1 would be more consistent with a 
10,000 yr old SNR. Our results, therefore, strengthen the case 
that the event recorded by Chinese astronomers was indeed a 
supernova and that RCW 86 is its remnant. 
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FIG. 4. — (Plate \)Chandra exposure corrected map of the northeastern part of RCW 86, using a square root brightness scaling. The red, green and blue channels 
correspond to the 0.5-1 keV 1-1.95 keV and 1.95-6.6 keV energy bands. 



